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ABSTRACT 

We present the results of a dedicated effort to measure the spin-down rate of the nearby isolated 
neutron star RX J0720.4— 3125. Comparing arrival times of the 8.39-sec pulsations for data from 
Chandra we derive an unambiguous timing solution for RX J0720.4— 3125 that is accurate to < 
0.1 cycles over > 5 years. Adding data from XMM and ROSAT, the final solution yields P = 
(6.98 ± 0.02) x 10~ 14 s s _1 ; for dipole spin-down, this implies a characteristic age of 2 Myr and a 
magnetic field strength of 2.4 x 10 13 G. The phase residuals are somewhat larger than those for 
purely regular spin-down, but do not show conclusive evidence for higher-order terms or a glitch. 
From our timing solution as well as recent X-ray spectroscopy, we concur with recent suggestions that 
RX J0720.4— 3125 is most likely an off-beam radio pulsar with a moderately high magnetic field. 

Subject headings: pulsars: individual (RX J0720.4— 3125) — stars: neutron — X-rays: stars 



1. INTRODUCTION 

One of the interesting results from ROSAT All-Sky 
survey ijVoges et al .111 9961 was the discovery of seven ob- 
jects that appear to be nearby, thermally-emitting neu- 
tron stars that have lit tle if any magnetospheric emis- 
sion (see HaberJ r2004l for a review). These objects, 
known most commonly as "isolated neutron stars," are 
distinguished by their long periods (> 3 s, when mea- 
sured), largely thermal spectra with cool temperatures 
(kT < 100 eV), faint optical counterparts (when de- 
tected), and lack of radio emission. 

It is not yet clear what sets the isolated neutron stars 
apart from the nearby, relatively young rotation-powered 
pulsars that also have cool thermal emission — sources 
like PSR B0656+14 and PSR B1055-52— which tend to 
have short (< 1 s) spin periods, ~ 10 12 -G magnetic fields, 
non-thermal (i.e. power-law) compon ents in their X- 
ray spectra, and rad io pulsations (e.g.. iPavlov & ZavlirJ 
120031 iKaplanl |2004) . The isolated neutron stars are 
known to have longer periods, but their spin-down rates 
(and hence magnetic fields ) are unknown, largely be- 
cause it has not yet been p ossible to determine a reliable, 
coheren t timing solution ijKaplan et alJ 120021 hereafter 
iPaner It fZane et al.ll20T)l . 

In this Letter we report a new analysis of the varia- 
tions of the 8.39-s period of the sec ond brightest sourc e 
of the group, RX J0720.4-3125 lIHaberl et alilT^ . 
We describe our analysis of Chandra data obtained 
specifically for timing purposes, as well as archival 
ROSAT, Chandra, and XMM data, in § In § we 
show that with the new data, we can avoid the pitfalls of 
the previous phase-coherent timing analyses and obtain 
a reliable timing solution. We discuss possible timing 
noise in § fl] and the implications of our result in § 
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2. OBSERVATIONS 

Our primary data were eight observations with 
the Advanced CC D Imaging Spectrometer (ACIS; 
iGarmire et, all2003T) aboard the Chandra X-ray Observa- 
torv < CXO : IWeisskonf et 111120001) . These were designed 
for timing accuracy, consisting of two sets of four expo- 
sures geometrically spaced over a period of about two 
weeks and separated by about half a year. We combined 
these with data from other Chandra observ ations, as well 
as fro m observation s with XMM-N ewton l).Iansen et alJ 
I2001|) and ROSAT ( Trumpet I\1)93Il A log of all obser- 
vations is given in Tabled 

For the Chandra data, we processed the level- 1 event 
lists to the level-2 stage following standard procedures 
and the latest calibration set (CALDB version 3.0.0). 
For the ACIS continuous-clocking data, this includes 
correcting the recorded event times for readout, dither, 
and spacecraft motion — corrections that used to re- 
quire additional steps ijZavlin et alJ I2000J) . We ex- 
tracted events within 1" of the source, and then ap- 
plied a clock correction of 284 .7 /zs to the arrival times 
ijDavis. Holmes. ~M vers 2003) ; the arrival times should 
now be accurate to < 6 (is. For the HRC-S/LETG data, 
we extracted zeroth-order events from a circle with ra- 
dius 10 pixels (1"3), and first-order events using the stan- 
dard LETG spectral extraction windows, but limited to 
10 < A < 65A. Finally, we used the axbary program 
to barycenter all of the events (using the optical po- 
sition: q. T ?nnn = 7 h 20 m 24!96, <5j 20 oo = -31°25'50"2; 
Kaplan et all 2003) 

For the XMM data, we used the standard procedures 
emchain and epchain (XMMSAS version 6.1.0) to repro- 
cess the observations. One additional step was necessary 
for the PN data set 622-U2, for which we found a small 
number of duplicate events (frames 963685-963719); we 
removed these before processing. Next, we extracted all 
single-pixel events within 37'.'5 of the source position, and 
used barycen to convert the arrival times to the solar- 
system barycenter 4 . 

4 Some portions of the 2000 and 2002 XMM /FN ob- 
servations were affected by a known processing problem 
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a PSPC: Position Sensitive Proportional Counter (Bricl & Pfcffcrmann 1995) 
aboard ROSAT. HRI: High-Resolution Imager (Zombcck et al. 1995) aboard 
ROSAT. HRC: High-Reso lution Camera for spectroscopy aboard Chandra 
(HRC-S; Kraft et al. 1997), used with the Low-Energy Tran smission Grat- 
ing. A CIS: Chandra's Advanced CCD Imaging Spectrometer (Garm ire et alj 
2003), used in continuous clocking mode. EPIC-pn: XMM's E uropean Pho- 
ton Imaging Camera with PN detectors jStriidcr et al. 2001), used in full- 
frame (ff) or small window (sw) mode, with thin, medium, or thick filter. 
EPIC-M0S1/2: European Photo n Imaging Cameras with MOS detectors 
aboard XMM (Turner et al. 2001), used in small-window mode with thin or 
no (open) filter. 

b Observation identifier (CXO, ROSAT) or revolution number and exposure 
identifier (XMM). 

c The TOA is defined as the time of maximum light closest to the middle of 
each observation, and is given with 1-a uncertainties. 



The r eduction of the ROSAT data followed that in 
iPaper % except that we properly corrected the event 
times to the Barycentric Dynamical Time (TDB) sys- 
tem instead of the Coordinated Universal Time (UTC) 
system returned by the FT00LS barycentering tasks (see 
iCronper et ^aU 120041) . We used the corrections supplied 



that rejected significant portions of the observations; see 
http : //xmm. vilspa. esa. es/sas/documentation/watchout/lost_e 
This should not introduce any systematic error, though it means 
that our TOAs for these observations are not as precise as possible 
with all events. However, since the present TOA uncertainties 
are smaller than the timing noise (§|1J, we decided not to try to 
remedy this problem. 



3. TIMING ANALYSIS 

Our goal was to use times-of-arrival (TOAs) to infer a 
phase-coherent timing solution involving the spin period 
and its derivative, where each cycle of the source was 
counted. To measure TOAs we needed an initial refer- 
ence period, something which we determine d using a Z\ 
test (Rayleigh statistic: lBuccheri et al.ll983|) on the com- 
bined 2004 January ACIS and 2004 February HRC data. 
We find P = 8.3911159(10) s (here and below, numbers 
in parentheses indicate twice the formal 1-a uncertain- 
slk&ntt- the last digit unless otherwis e indicat ed), which 
is consistent with our earlier value l|Paner ff l but much 
more accurate because we could coherently connect ob- 
servations over a much longer (52 day) time span. 

Using this period, we constructed binned light curves 
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MJD (days) x10 < 

Fig. 1.— Phase residuals for RX J0720.4-3125. The top panel 
shows the residuals for each TOA compared to a linear (y = 0) 
model. The solid curve gives the best-fit quadratic (u ^ 0, v = 0) 
ephemeris for all data (Tab. [2] with a 0.27-s systematic uncertainty 
added in quadrature). The vertical dashed line indicates the ref- 
erence time to- The bottom panel shows the residuals relative to 
the quadratic model. We also show a best-fit cubic model {y ^ 0; 
dotted line) and a model that includes a glitch in frequency near 
MJD 52800 (A/ « 1 x 10~ 9 ; dashed curve). 



(with 16 phase bins) and determined the TOAs by fitting 
a sinusoid (approp riate for the sinus oidal pulsations of 
RX J0720.4-3125; Hab erl et al1ll997|) : the uncertainties 
were calculated from the uncertainties in the phases of 
the fitted sinusoids (we verified that we obtain consistent 
results if we change the binning or measure TOAs using 
cross-correlation instead). Here, the TOA is defined as 
the time of maximum light closest to the middle of the 
observation, a choice which minimizes co-variance with 
small changes in period. We present TOAs for all of the 
data in Tabled 

We then determined a timing solution for only the 
Chandra data using an iterative procedure. We first used 
the reference period to determine cycle counts for the five 
above-mentioned observations, as well as the next obser- 
vation closest in time. We fit these cycle counts to 

4>(t) = 0o + v{t - 1 ) + 2^(* - + g v(t - 1 ) 3 • • • , (1) 

where <po is the cycle count plus phase at reference time 
to, v is the spin- frequency, v is its derivative, v is the 
second derivative. We then iterated, using the improved 
solution to determine the cycle count for the next ob- 
servation, etc. We started with v = v = 0, but left v 
free once that significantly improved the fit; v was not 
required (cf., §0] below). 

The final ephemeris listed in Table [21 has small, < 
0.1 cycle, residuals, and fits the Chandra data well: 
xl = x 2 /Ndof — 1-06 (with A do f = 13 degrees of free- 
dom). To test the uniqueness of our ephemeris, we tried 
changing the cycle counts (adding or subtracting one or 
more cycles) at the least unambiguous points, but found 
that the resulting solutions were very poor (e.g., alter- 
ing the cycles between the 2001 ACIS and 2000 HRC 
observations gave xl = 89.26). 

To improve and extend the ephemeris, we added the 
XMM and ROSAT data into the solution (see Table EJ. 
As for the Chandra data alone, the cycle counts are 
unambiguous and, as can be seen in Fig. ^ the resid- 



TABLE 2 

X-ray Timing Ephemerides for RX J0720. 4-3125 



Quantity a CXO only All data b 



Dates (MJD) 51575-53394 49257-53394 

to (MJD).... 53010.2635605(18) 53010.2635637(24) 

v (Hz) 0.11917366926(46) 0.11917366908(38) 

v (Hz s" 1 ) . . -9.97(11) x KT 16 -9.918(30) x 10~ lf 

TOA rms (s) 0.18 0.31 

X 2 /DOF .... 13.8/13=1.06 30.7/31=0.99 

P (s) 8.391115305(32) 8.391115532(26) 

Piss- 1 ).... 7.019(80) x 10" 14 6.983(22) x 10" 14 

B (ergs" 1 ).. 4.7 x 10 30 4.7 x 10 30 

B dip (G).... 2.5 x 10 13 2.4 x 10 13 

7-char (yr) .... 1.9 x 10 6 1.9 x 10 6 



Note. — Uncertainties quoted are twice the formal 1-cr 



uncertainties in the fit. 

a _E = 10 45 /4547r 2 ^!> is the spin-down luminosity (with 
/ = 10 45 745 g cm 2 the moment of inertia); £?dip = 3.2 x 

10 19 Vpp is the magnetic field inferred assuming spin-down 
by dipole radiation; T cllar = P/2P is the characteristic age, 
assuming an initial spin period Po <S P n ow 

b A 0.27-s systematic uncertainty has been added in 
quadrature to these data to account for timing noise. See 

□ 



uals remain below 0.1 cycles, lending additional cre- 
dence to our results. For verification, we also exam- 
ined the fit using just the Chandra and XMM data, 
which avoids the large gaps between the ROSAT obser- 
vations. We found the same cycle counts, but a statisti- 
cally significant difference in the inferred parameters [v — 
0.11917366887(12) Hz and v = -9.69(4) x 10~ 16 Hz s _1 , 
while the fit to all data gave v = 0.1191736700(2) and 
v = -9.915(14) x 10" 16 Hz s" 1 ]. 

While we arc confident in our fit in general, the above 
discrepancy is puzzling. Furthermore, the xl values for 
the fits including the XMM data are poor, with xl = 
6.16 for Chandra+XMM (A dof = 28; rms = 0.34 s) and 
xl = 10.25 for all of the data (N dof = 31; rms = 0.36 s), 
while the solution for Chandra+ROSAT, although some- 
what ambiguous, is tolerable (xl = 1.40) and similar to 
the Chandra-only solution. The x 2 values for the full fits 
are dominated by the contribution of the XMM data, 
which have very small formal uncertainties (but are not 
always entirely consistent from one instrument to an- 
other for the same observation; see Tab.^l. The XMM 
data also cause the values of i> to differ by 3% and the 
Chandra+XMM ephemeris to be a poor match to the 
ROSAT points. We discuss these deviations in more de- 
tail below; here, we note that for the estimates of the val- 
ues and uncertainties in Table we ensured xl — 1 by 
adding in quadrature an additional uncertainty of 0.27 s 
to all TOAs. We stress, however, that our overall solu- 
tion is robust, and the inferred spin-down rate should be 
reliable at the < 3% level. 

4. DEVIATIONS FROM REGULAR SPIN DOWN? 

There could be several reasons for the relatively poor 
fit of the timing solutions to our full set of arrival times. 
First, RX J0720.4— 3125 may show some rotational insta- 
bilities or "timing noise," as seen in radio pulsars. We 
can estimate the magnitude of the long-term timing noise 
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by fitting the phase residuals with a third-order [p ^ 0) 
solution like that shown in Fig. ^ This does in fact im- 
prove the fit — giving xl = 6.02 for N^of = 30 with an 
rms of 0.32 s — though it is still formally unacceptable. 
We find v 2.4(4) x 10~ 25 Hz s~ 2 (this also changes 
v by 3% from the value in Tab. giving a timing 
noise measurement of A 8 = log 10 (|i>| (10 s ) 3 /6i^) = —0.5. 
This value is on the high side of, but not outside the 
range expected fro m relations between P an d Ag found 
for rad io pulsars l|A rzoumanian et a lJ 119941) and mag- 
netars l|Woods et all 120001: iGavriil fc Kaspil l2002j) . We 
do not believe it is likely that the third-order solution 
represents a real long-term change in i>, since the value 
of v changes significantly if one uses just the Chan- 
dra, Chandra+XMM, or all of the data 5 [8(6), 18(4), 
or 2.4(4) x 10~ 25 Hz s~ 2 , with xl = 0.64, 2.81, or 6.02 
respectively] . 

A second possibility is a sudden change in rota- 
tion — a glitch, as proposed (and largely rejected) by 
Ide Vries et al.l ()2004|) to account for variations in spec- 
tral shape (differences of 10% in the inferred temper- 
ature) and pulse shape (changes of 50-100% in pulsed 
fraction) observed between 2000 and 2003. The system- 
atic increase in the residuals after MJD 53000 in Fig- 
ure n might indeed indicate a glitch. We tried fitting a 
simple glitch model, in which we assumed that only the 
frequency changed, that the recovery time was longer 
than the span of our observations, and that the glitch 
occurred on 2003 July 1 (MJD 52821), in between the 
two XMM observations that sh owed the largest spec- 
tral change l)de Vries et alJ l2004'). We find a reasonable 
fit for a glitch with A/ = 1.3 x 10 -9 Hz and a recov- 
ery time > 3 yr, giving xl — 7.7 (see Fig. [TJ. This 
would be a small glitch, with A/// = 1.1 x 10 -8 com- 
pared to values of 10 -9 to 10~ 6 for radio pulsars (with 
smaller values more typical for larger magnetic fields; 
iLvne. Shemar. fc Smithll2000|) . It also implies that ener- 
getically, the putative glitch would be insig nificant: the 
chang e in kinetic energy of ~ 10 36 ergs l|van Riper et al.l 
1991) would only be noticeable if dissipated in < 1 day 
given the bolometric luminosity of RX J0720.4— 3125 of 
2 x 10 32 d 2 00 ergs s _1 (whe re d = 300d3on pc i s the dis 



in pc i s 

tance to RX J0720.4-3125: lKaplan et al.l2003j) . In prin- 
ciple, however, it could still have altered the light curve 
and spectrum of RX J0720.4— 3125 through realignment 
of the magnetic field relative to the spin axis. 

Finally, a more mundane explanation for the rela- 
tively poor fit is that the data are from different instru- 
ments with different energy responses (even among a sin- 
gle instrument aboard XMM, the changing filters alter 
the response). The pulse profile of RX J0720.4-3125 
is known to depend on en ergy l)Cropper et alJ 1200 H 
IPaper ft IHaberl et all I2004|) and to change over time 
<jde\^eset^dTl2"o(Mjr ^ While the changes to the shape 
are small, some systematic offsets are expected between 
the pulse profiles as measured by different instruments or 
at different times. We hope to investigate this in more 
detail in the near future. 

At present, we cannot distinguish between the various 

5 From Fig.0 it may appear that one could obtain a better fit by 
reducing the cycle count for the 1998 HRI point by one; however, 
doing that, the other ROSAT points cannot be reproduced any 
more. 



possibilities for the relatively poor fits. The predicted 
future behavior is different, however, and thus further 
observations of RX J0720.4— 3125 (some of which are 
in progress) should be able to distinguish between these 
models. 

5. DISCUSSION & CONCLUSIONS 
From our timing solution, we infer a spin-down rate 
P = 6.98(2) x I P" 14 s s" 1 . Thi s is consistent wit h the 
limits derived in IPaper J and by IZane eTaf] 112001, but 
inconsistent with the tentative solution of lCropper et all 
i200l . who found P = (1.4 ± 0.6) x 10" 13 s s" 1 at 99% 
confidence (but who noted that elements of their solution 
were inconsistent with each other and that their analy- 
sis was subject to confusing aliases). In Table [21 we list 
derived parameters — rotational energy loss rate, mag- 
netic field, and characteristic age (the latter two under 
the assumption of magnetic dipole spin-down). 

The values of P and P place RX J0720.4-3125 well 
above above mos t of the radio-pulsar "death-lines" pro- 
posed so far (e.g.. lYoung. Manchester, fc Johnston! 1999J) 
and in a region popu lated b y radio pulsars in P-P 
diagrams like that in IPaper il (its parameters are ap- 
proximately between those of PSR J1830— 1135, with 
P = 6.2 s, P = 5 x 10~ 14 s s-\ and PSR J1847-0130, 
with P = 6.7 s, P = 1.3 x 10~ 12 s s" 1 ). Hence, 
RX J0720. 4—3125 may well be a radio pulsar itself, but 
one whose radio beam(s) do not intersect our line of sight. 
Its inferred magnetic field, B — 2.4 x 10 13 G, is not excep- 
tiona l; the Parkes Multi-beam Survey ({Manchester et all 
|2001[) in particular has discovered a fair number of ra- 
dio pulsars with B > 10 13 G (e.g.. JCamilo et alJl2000l 
IMorris et alJ2002tlMcLaughlin et alJ2003fl . and it is now 
clear the distributio n of magnetic fields is flatter than 
previously assumed l)Vranesevic et alj|2004|) . 

With E = 4.7 x 10 30 ergs s~\ RX J0720.4-3125 is 
not expected to have muc h non-thermal X-ra y emis - 
sion: from the relation of iBecker fc Trumpet! l)1997|) . 
one estimates Lx,non-th ~ 10~ 3 E ~ 5 x 10 27 ergs s _1 , 
much smaller than the thermal emission, Lx, therm — 
2 x 10 32 d3oo ergss -1 . This is consistent with lim- 



its from Chandra and XMM llPaerels et alJ 1200 H 
IPavlov. Zavlin. fc Sanwalll2002t IRaplan et alJl2003ft . 



What is somewhat puzzling is the inferred age of 2 Myr. 
Tracing RX J0720.4— 3125 back to OB associations where 
it might have been born put it close to the Trumpler 10 
assoc i ation ~ 0.7 Myr ago l)Motch. Zavlin. fc Haberl 
120031 lKaplanl[2004|) . Similarly, based on its estimated 
temperature and luminosity, most standard cooling mod- 
els (modified URCA for 1.4 M w neutron stars) put 
RX J0720.4-3125 at < 1 M yr lIHevl fc HernauistifToM 
IPaper It ICropper et alJl200^l . 

It is of course possible that RX J0720.4— 3125 was born 
with a long period and/or had significantly non-dipole 
braking, such that the spin-down age is not a good es- 
timate of its true age. However, no case with as long a 
birth period as would be required f or RX J0720.4— 3125 
is known among r adio pulsars (cf. IKramer et alJ 120031 
IGavriil et aDl200l . 

Another possible explanation is that RX J0720.4— 3125 
was ejected from a binary system with a massive com- 
panion ~ 0.7 Myr ago, either when the companion un- 
derwent a supernova or during a binary exchange inter- 
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action. In this relatively long period is expected: 

if the neutron star accreted matter from its compan- 
ion, its spin period would have tended toward the equi- 
librium period, F cq ps 5 s (B/10 13 G) 6 / 7 (M/A/ E dd)~ 3/7 
(where M is the accretion rate and MEdd is the Edding- 
ton rate). A relatively short cooling age would also be 
consistent with this model: the accretion and accompa- 
nying steady hydrogen burning could reheat the neutron 
star (or keep it hot). Of course, it remains to be seen that 
a suitable evolutionary scenario can be found. In any 
case, the prediction for the model with two supernovae 
is that there may well be another < 1 Myr old neutron 
star whose proper motio n traces back to the sam e origin 
as RX J0720.4-3125 (cf. IVlemmings et alll200l|) . 

Finally, we can compare the magnetic field strength of 
2.4 x 10 13 G with what is inferred from the broad absorp- 
tion feature in the spectrum (observed to be at 0.3 keV, 
which corresponds to 0.3 9 keV at the surface for a grav- 
itatio nal redshift of 0.3; lHaberl et al.l 120041 iVink et all 
I2004J) . If due to a proton cyclotron line, one infers 



B~6x 10 13 G ijHaberl et aLl UOQ^. which is substan- 
tially larger than inferred from the spin-down. This may 
simply reflect the inadequacy of the dipole spin-down 
model, or the presence of higher order multipoles. On the 
other hand, based on a c omparison with other sources, 
Ivan Kerkwiik et all l)2004l) suggested that the absorption 
feature was due to the transition from the ground state to 
the second excited tightly bound state of neutral hydro- 
gen, which would require B ~ 2 x 10 13 G and matches the 
spin-down value nicely. If that is the case, higher signal- 
to-noise spectra should reveal the transition to the first 
excited state at ~0.15 keV. 
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